Journal of the Physical Society of Japan 



FULL PAPERS 



o 

(N 
Or 
C/3 



a,: 

-I— > 



> 
O 

o 

m 



X 



Aperture 



Scalable Neutral Atom Quantum Computer with Interaction on 
Demand: Proposal for Selective Application of Two-Qubit Gate 
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We propose a scalable neutral atom quantum computer with an on-demand interaction through a selective 
two-qubit gate operation. Atoms are trapped by a lattice of near field Fresnel diffraction lights so that 
each trap captures a single atom. One-qubit gate operation is implemented by a gate control laser beam 
which is applied to an individual atom. Two-qubit gate operation between an arbitrary pair of atoms is 
implemented by sending these atoms to a state-dependent optical lattice and making them collide so that a 
particular two-qubit state acquires a dynamical phase. We give numerical evaluations corresponding to these 
processes, from which we estimate the upper bound of a two-qubit gate operation time and corresponding 
gate fidelity. Our proposal is feasible within currently available technology developed in cold atom gas, 
MEMS, nanolithography, and various areas in optics. 
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1. Introduction 

Quantum computation employs a quantum mechan- 
ical system as a computational resource. 1-3 ) It is ex- 
pected that a large scale quantum computer outperforms 
its classical counterpart exponentially by making use of 
superposition states and entangled states. In spite of 
many proposals for potentially scalable quantum com- 
puters, most physical realizations so far accommodate 
qubits on the order of ten, at most. One of the obsta- 
cles against scalability in the previous proposals is the 
absence of controllable interaction between an arbitrary 
pair of qubits. Small scale quantum computers are al- 
ready demonstrated experimentally in several physical 
systems, such as NMR, trapped ions, photons, and neu- 
tral atoms. These physical systems have merits and de- 
merits and no single physical system to date satisfies all 
the DiVincenzo criteria, the necessary conditions for a 
physical system to be a candidate of a working quantum 
computer. 4 -' 

A neutral atom quantum computer is one of the most 
promising candidates of a scalable quantum computer. 
Neutral atoms are believed to be robust against decoher- 
ence, another obstacle against a physical realization of a 
working quantum computer, since they are coupled very 
weakly to the environment. A one-qubit gate operation 
has been already demonstrated with a two-photon Ra- 
man transition, 5 ' 6 ) which is a well-established technique 
today. A two-qubit gate operation makes use of a cut- 
ting edge technique involving ingenious manipulation of 
hyperfine-state-dependent optical lattice potentials. 7 ~ 10 ) 
A drawback of this two-qubit gate implementation is that 
the gate acts on all the nearest neighbor pairs in the op- 
tical lattice simultaneously. Hence a selective gate oper- 
ation is yet to be realized, although some ideas utilizing 
optical tweezers has been proposed. n ) 
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Fig. 1. (Color online) (a) An atom trapped in an NFFD light. 
Red arrow (solid line) shows the trap light while the green arrow 
(broken line) shows the laser light required for one-qubit gate op- 
erations, (b) Example of an array of such traps. The trap laser and 
the gate control laser are explicitly shown only for the first row. (c) 
A one-dimensional optical lattice is superimposed near the array 
of the NFFD traps. The minima of the optical lattice are situated 
so as to adjust the period of the array of the NFFD traps. 



It is the purpose of the present paper to propose a new 
design of a neutral atom quantum computer with an on- 
demand interaction, which is potentially scalable up to 
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nology. 12, 13 ) The schematic illustration of our proposal 
is shown in Fig. 1. We punch an array of apertures in a 
thin substrate made of, e.g., silicon. An optical fiber is 
attached to each aperture, which lets in two laser beams, 
one for trapping and the other for one-qubit gate control. 
Atoms are trapped by a near field Fresnel diffraction light 
(NFFD), 14 ) which is produced by letting the trap laser 
light pass through an aperture with a diameter compara- 
ble to the wavelength of the laser, so that each trap holds 
a single atom. Turning on and off the trap laser individ- 
ually is required to implement a selective two-qubit gate 
as will be shown later. Each lattice site is also equipped 
with its own gate control laser beam so that individ- 
ual singlc-qubit gate operations can be applied on many 
qubits individually and simultaneously. A controlled two- 
qubit gate operation between an arbitrary pair of qubits 
is made possible by selectively transferring two atoms in 
a superimposed one-dimensional optical lattice and by 
employing the method demonstrated by Mandel et al. 9 ^ 
Here, we give detailed analysis of the selective two-qubit 
gate implementation according to our proposal. 

Technology required for physical implementation of 
the proposed quantum computer has already been devel- 
oped in the areas of cold atom trap, scanning near field 
optical microscope (SNOM), nanolithography, and quan- 
tum optics, among others, and we believe there should 
be no obstacles against its physical realization. 

We first briefly review the design of a neutral atom 
quantum computer demonstrated by Mandel et al. 9 ^ in 
Sec. 2. We point out problems inherent in their design. 
Individual atom trap by making use of NFFD light is 
introduced in Sec. 3. Section 4 is devoted to one- and 
two-qubit gate implementations. We separate the two- 
qubit gate operation into several steps and evaluate the 
execution time and the corresponding fidelity of each step 
in Sec. 5 and Appendix A. Conclusion is given in Sec. 6. 

2. Brief Review of Neutral Atom Quantum 
Computer in Optical Lattice 

In this section, we look back the experiments by Man- 
del et alJ~ 9 ^ to give a basic background on how to trap 
and manipulate atoms using laser beams and how to im- 
plement gate operation. We subsequently point out prob- 
lems inherent in their design. 

Suppose an atom is put in a laser beam with an os- 
cillating electric field E(x,t) = Re(E (x)e~' lULt ), where 
ujl is the laser frequency. It is assumed that w_l is close 
to some transition frequency ujq = E e — E g between two 
states \g) and |e) of the atom. The interaction between 
the electric field and the dipole moment of the atom in- 
troduces an interaction Hamiltonian 

H i = ~(E -dC)(e- iULt + e iui - t ), (1) 

where d is the dipole moment operator of the atom. This 
interaction introduces an effective potential of the form 

m _ H^pl m 

for an atom in the ground state, which is called the 
AC Stark shift. Here il eq = (e\d\g) ■ E (x)/h, while 



A eg = ujl— ujq is the detuning. We have ignored the small 
natural line width of the excited state. In case A eg > 
(blue-detuned laser), V(x) is positive and a region with 
large V(x) works as a repulsive potential. In contrast, 
if A eg < (red-detuned laser), V(x) is negative and a 
region with large | V(as)| works as an attractive poten- 
tial. We use mainly a red-detuned laser in the following 
proposals. 

It is possible to confine neutral atoms by introducing 
a pair of counterpropagating laser beams with the same 
frequency, amplitude and polarization along the x-axis. 
These beams produce a standing wave potential, called 
an "optical lattice" of the form 

Vol(x) = -Vocos 2 (k OL x)e^y 2+z2 y w \ (3) 

where fcoL is the wave number of the laser and w is the 
waist size, which is on the order of the laser wavelength. 
It is possible to construct a three-dimensional lattice by 
adding two sets of counterpropagating beams along the 
y- and the z-axes. Note that the lattice constant is always 
Aol/2, where Aol is the wavelength of the laser. 

Mandel et at trapped atoms in such an optical lat- 
tice. 9 ) They used the Rabi oscillation to implement one- 
qubit gates, in which a microwave (MW) field was ap- 
plied. They have chosen qubit basis vectors |0) = \F — 
l,m_p = —1) and |1) = \F = 2,rriF = —2) to be in har- 
mony with the one-qubit gate operation making use of 
the Rabi oscillation. 

For two-qubit gate operations, they introduced time- 
dependent polarization in the counterpropagating laser 
beams as 

E+(x) = e tk °^ x (zcos9 + ysm9), 

(4) 

E_(x) = e' lkoLX (z cos (9 -y sin (9). 

These counterpropagating laser beams produce an opti- 
cal potential of the form 

E + (x) + E-(x) oc er + cos(fcoLa; — 9) + o~ cos(/col£ + 9), 

(5) 

where <r + (a ) denotes counterclockwise (clockwise) cir- 
cular polarization. Note that the first component (oc a + ) 
moves along the x-axis as 9 is increased, while the sec- 
ond component (oc cr~) moves along the —x-axis un- 
der this change. 8 ^ The component a + introduces tran- 
sitions between fine structures; nSi/ 2 (mj = —1/2) — > 
nPi/2 (mj = 1/2), nSi/2 {mj = -1/2) -> nP 3/2 (mj = 
1/2), and nSi/a (mj = 1/2) -> nP 3/2 (m, ; = 3/2). 
The transitions from nS\/ 2 to nP 3 / 2 are red-detuned, 
while the transition from nSi/ 2 to nP 1 / 2 is blue-detuned 
if ul is chosen between the transition frequencies of 
nSi/2 (mj = —1/2) — > nP 3 / 2 (mj = 1/2) and 
nS 1/2 (mj = -1/2) -> nP 1/2 (mj = 1/2). Then by 
adjusting ujl properly, it is possible to cancel the at- 
tractive potential and the repulsive potential associated 
with these transitions. The net contribution of the a + 
laser beam in this case is an attractive potential for an 
atom in the state nSi/ 2 (mj = 1/2). We denote this 
potential as V + (x) oc cos 2 (fcoL^ — 9) in the following. 
Similarly, the <r~ component introduces a net attractive 

potential V_(x) oc cos 2 (koi J x+9), through the transition 

„ c . /™ 1 /o\ k „ D /™ 7 _ q /o\ ™ 
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in the state nSi/ 2 {m,j = —1/2). 

Mandel et a/. 9 ' took advantage of these state- 
dependent potentials to implement a two-qubit gate. 
By applying the Walsh-Hadamard gate on |0)j|0)j+i, 
one generates a tensor product state (|0) + |l))j(|0) + 
|l))j + i/2. The potentials acting on the states |0) and |1) 
are evaluated as 

V lQ) (x) = lv + (x) + \v^x) (g) 

V\ 1} (x) = V-(x). 

By decreasing the phase 9, the state |0) moves with dom- 
inating V+ toward -l-cc-direction, while |1) moves with V— 
toward — x-direction. Thus it is possible to make |0)j and 
|l)i+i collide between the two lattice points. If they are 
kept in the common potential well during ihoid, the sub- 
space |Q)i|l)j+i obtains a dynamical phase e _lC/int * hold , 7 ) 
where Ui n t is the on-site repulsive potential energy. After 
these two components spend thoid in the potential well, 
they are brought back to the initial lattice point by re- 
versing 9, which results in the two-qubit gate operation 

|0)i|l) i+1 V|0) i |0) i+1 + e- i ^«|0) i |l) i+1 

+|l) j |0) i+1 + (7) 

It should be noted, however, that the state dependent 
potentials act on all the pairs of the atoms and selective 
operation of the two-qubit gate on a particular pair is im- 
possible. Their proposal may be applicable to generate a 
highly entangled state for a cluster state quantum com- 
puting, although it is not applicable for a circuit model 
quantum computation. 

In the following sections, we propose implementations 
of a neutral atom quantum computer, which overcome 
these difficulties. 

3. Near Field Fresnel Diffraction Trap 

Suppose there is an array of apertures in a thin sub- 
strate. We consider a square lattice of apertures for def- 
initeness. The substrate is made of a thin silicon, for 
example. Silicon has a very good thermal conductivity, 
comparable to a metal, and local heating due to a laser is 
expected to diffuse quickly. It is possible to maintain the 
substrate in a low temperature if the perimeter of the 
substrate is in contact with a refrigerator, which sup- 
presses thermal excitation of an atom. 

Bandi et al. 1A ^ proposed to trap an atom with a mi- 
crotrap employing NFFD light. Atoms are trapped in an 
array of NFFD traps, each of which traps a single atom. 
An NFFD light is produced if a plane wave with the 
wave length \p is incident to a screen with an aperture 
of the radius a > Xf. The ratio Np — a/Xp > 1 is called 
the Fresnel number. The trap potential is evaluated by 
applying the Rayleigh-Sommerfeld formula as 

U F (x) = -t/o^S^, (8) 




Fig. 2. Potential profiles of the NFFD trap for (a) a/A F = 1.5, 
and (b) a/Xp = 2.8, where a is the aperture radius and Ap the 
wavelength of the trap laser beam. The left panel is the contour 
plot while the right panel shows the potential profile along the 
z-axis with x = y = 0. The length and the potential energy are 
normalized by \p and Uq, respectively. 



where 

S(x) = ^JJ^-^(l- ik ?) dx ' d v' ( 9 ) 

and 

with the distance r between x = (x,y,z) and x' = 
(x 1 , y' , 0) in the plane. Here, T e is half of the spontaneous 
decay rate, E$ the amplitude of the incoming plane wave, 
kp = 2-k/Xp its wave number, the detuning A eg is neg- 
ative (red-detuned), and the plane wave is incident to a 
screen from the —z-axis. The integral is over the aperture 
region (x 12 + y' 2 < a 2 ). 

The aperture size is on the order of the wave length 
so that the distance between the minimum of the po- 
tential energy and the substrate is on the order of the 
wavelength. Figure 2 shows the trap potential profiles 
for two aperture radii, which is obtained by evaluating 
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Fig. 3. z-coordinate of the potential minimum (z m ) as a func- 
tion of the aperture size a/Ap. 



Eq. (8) numerically. It is clear from Fig. 2. that the dis- 
tance between the atom and the substrate is controllable 
by manipulating the aperture size a. Figure 3 shows the 
z-coordinate of the potential minimum (z m ) as a func- 
tion of the aperture size a. We take advantage of this fact 
in our implementation of a two-qubit gate later. 

The aperture size may be controlled by employing a 
small scale liquid crystal technology called Spatial Light 
Modulator (SLM) with Liquid Crystal on Silicon (LCOS) 
technology. 15 ' 16 - 1 An array of shutters made of liquid 
crystal may be fabricated on a silicon wafer with LCOS- 
SLM. The current operation time of the SLM is reported 
to be sub msec, which is comparable to the operation 
times of the other steps in our gate implementation 
discussed below. Alternatively, well-established Micro- 
Electro-Mechanical Systems (MEMS) technology may be 
employed to control the aperture radius. An MEMS shut- 
ter for a display is already in a mass production stage 17 -* 
and it may be applied to demonstrate our proposal. The 
on (off) switching time of the shutter is 54 (36) /is 17 ' 
and we expect even faster switching time by optimizing 
its structure for our purpose. 18 ' 

It is clear by construction that we have freedom in the 
choice of the lattice constant and the lattice shape. We 
may even arrange the traps in an irregular form if nec- 
essary. This is a remarkable new feature, which a con- 
ventional two-dimensional optical lattice cannot provide 
with; an optical lattice is always regular and its lattice 
constant is on the order of the wavelength of the laser 
beams, where the lattice constant can be tuned to some 
extent by changing the angle of two lasers. 19 ' 

4. Quantum Gate Implementations 

For a quantum system to be a candidate of a general 
purpose quantum computer, it is necessary for the sys- 
tem to be able to implement a universal set of quantum 
gates, that is, the set of one-qubit gates and almost any 
two-qubit gate. 20 ' Now we explain how these gates are 
realized in our proposal. 

In the following we take two hyperfine states to repre- 
sent two-qubit basis vectors, which we denote |0) and |1). 
The hyperfine states must be consistent with two-qubit 



1 , txi p = 1) and |1) = \F = 2,mp = 1) throughout this 
paper. 

4-1 One-Qubit Gates 

Let us assume, for definiteness, the NFFD trap poten- 
tial is strong enough so that there is only one atom in 
each trap. Along with the NFFD trap light, there is a 
gate control laser beam propagating through the fiber, 
which is depicted in a broken line in Fig. 1. 

One-qubit gate is implemented by making use of the 
two-photon Raman transition, 3 ' which is already demon- 
strated. 5 ' 6 ' Let Eq and E\ be the energy eigenvalues of 
the states |0) and |1), respectively, and let E e be the 
energy eigenvalue of an auxiliary excited state |e) nec- 
essary for the Raman transition. Suppose that a laser 
beam with the frequency uj^ has been applied to the 
atom. Let A = u L — (E e — E )/h be the detuning 
and fij be the Rabi oscillation frequency between the 
state \i) (i — 0,1) and |e). Then, under the assump- 
tions |A| 3> (Ei — Eo)/h, fi 2 /|A|, we obtain the effective 
Hamiltonian 

Hi = -e<r z - (11) 

where 

e = Ei-E + ^°. 

Note that Hi generates all the elements of SU(2) since 
there are two su(2) generators cr XtZ in the Hamiltonian 
and their coefficients are controllable. 

4-2 Two-Qubit Gates 

As illustrated in Fig. 4, our proposal for selective two- 
qubit gate operation consists of 7 steps: 

1 We choose two trapped atoms to be operated by the 
gate. Each atom is made into a superposition state 
(|0) + \l))/y/2 by a gate control laser beam attached 
to each aperture. Then, these atoms are transferred 
from the vicinity of the substrate to space coordi- 
nates, where a one-dimensional optical lattice will be 
introduced later, by enlarging the aperture radius. 
The minima of the optical lattice will be situated at 
these space coordinates. 

2 While the atoms stay at the points, the NFFD lasers 
of the two atoms are gradually turned off while a 
pair of counterpropagating laser beams to form an 
optical lattice is gradually turned on. 

3 Now the two atoms are trapped in the optical lattice. 
As described in Sec. 2, the polarizations of the pair 
of the counterpropagating laser beams are rotated 
in opposite directions so that the qubit state |0) is 
transferred toward positive x-direction, while |1) is 
transferred toward negative direction. Therefore, it 
is possible to collide |0) of one atom and 1 1) of the 
other atom. 

4 Now, component |0) of one atom and |1) of the 
other are trapped in the same potential well in 
the optical lattice. They interact with each other 
for a duration tu^A so that this particular state 
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|0)|1) acquires an extra phase J/mtihoid compared 
to other components |0)|0), |1)|0), and |1)|1). Here 
Uint — (4:Trh 2 a s /m) J \ip\ 4 dx is the on-site interac- 
tion energy with the atomic mass m and the s-wave 
scattering length a s . When thoid satisfies the condi- 
tion ?7i n tihoki = 7r, we obtain a nonlocal two-qubit 
gate |00)(11| -|01)(01| + |10)(10| + 111) (11|. Applica- 
tion of further one-qubit gates implements a CNOT 
gate or a CZ gate. 

3' After acquiring the phase, two atoms are separated 
along the optical lattice by an inverse process of Step 
3. 

2' Atoms are transferred from the optical lattice to the 
NFFD trap by an inverse process of Step 2. 

1' Aperture radii of NFFD lasers are reduced so that 
the atoms are transferred back to the vicinity of the 
substrate by an inverse process of Step 1. Now the 
gate operation is completed. 

5. Execution Time of Two-Qubit Gate 

In this section, we analyze all steps of selective two- 
qubit gate operation in detail one by one to estimate 
the execution time. While shorter execution time is pre- 
ferred, the fidelity should be kept close to 1 during the 
gate operation. The fidelity is mostly bounded by adi- 
abaticity requirement. Here, we give an estimate of the 
execution time, which is limited by the adiabatic con- 
dition. Detailed numerical simulations, using the time- 
dependent Schrodinger equation with experimentally re- 
alistic parameters, for each step will be described in Ap- 
pendix A. 



□Step 1 ' 



Fig. 4. (Color online) Schematic diagram of a two-qubit gate 
operation. In the top panel, the left and the right atoms are acted 
by the Walsh-Hadamard gate in advance so that each of them is in 
the superposition of |0) (red semicircle) and ]1) (green semicircle). 
Here, the black contours show the individual NFFD potential and 
bold black segments show the screen. Step 1: By enlarging the 
aperture size, the minima of the selected traps are moved away 
from the substrate. Step 2: A one-dimensional optical lattice along 
the x-axis, depicted by the thin red contours, is superimposed with 
turning off adiabatically the selected NFFD traps. Two atoms in 
the superposition of |0) and |1) are left in the state-dependent 
optical lattice. Step 3: The polarizations of the counterpropagating 
laser beams are tilted by ±0 so that the |0) and |1) components 
of two atoms move in the opposite directions. Step 4: The collision 
of |0) of one atom and |1) of the other gives a dynamical phase. 
Steps 3', 2', and 1' are the reverse process of steps 3, 2, and 1, 
respectively. 



5.1 Step 1 

Suppose that an atom is transferred away from the 
vicinity of the substrate (from z m = z ini to zg n ) by in- 
creasing the aperture radius a, e.g. from (a) to (b) of 
Fig. 2. The distance Zf\ D — z ini may be a few times of 
Xf- As seen from Fig. 2, the NFFD trap potential in this 
parameter range is shallow along the z-axis and steep in 
the xy-plane, which means that the curvature along the 
z-axis dominates the adiabaticity condition. The angular 
frequency along the z-axis around the minimum z — z m 
is 



1 d 2 U ¥ 



m dz 2 



\ 



2a- 



d 2 U* 



dz 2 



(12) 

TF 



where we put tildes for dimensionless variables scaled by 
the wavelength Ap and the potential amplitude Uq of 
the NFFD trap Eq. (8); rp = H/Uo is the characteristic 
time scale for NFFD trap and a = h 2 /2mf/oA F . The 
dimensionless curvature E/ F '(z m ) = d 2 Up/ dz 2 \z=s m as a 
function of the aperture size a is shown in Fig. 5. 

Since lu z is a monotonically decreasing function of a, 
the characteristic time t£j for an adiabatic change can 
be estimated from the condition 
/ \ 2 

1 / -2-fin -2-ini \ 



-m 



r (1) 

'ad 



hw z {z ri 



ZfinJ 



(13) 
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Fig. 5. (Color online) Curvature of the potential along the z-axis 
at the potential minimum as a function of the aperture radius. The 
length and the energy are normalized by and Uq, respectively. 



which gives 



r (1) 

'ad 



Zfin 



_ r „ - TT Tf (14) 

Here we have taken the "worst" curvature at the final 
value of a, i.e., z m = zr u . For the typical parameter val- 
ues for a realistic NFFD trap (see AppendixA.l) we can 
see a C 1, because the depth of the NFFD trap is very 
deep. We expect, from this estimation, that the trans- 
fer of an atom from the initial position to the optical 
lattice potential minimum will be bounded in time by 



lO^-Vp ~ 10 2 "Vs 



5.2 Step 2 

Now that the atoms are transferred to the right space 
points, their wavefunctions must be transformed from 
the ground states of the NFFD traps to those of the 
optical lattice adiabatically. It is important to note that 
the centers of the NFFD traps and those of the optical 
lattice potential must be sitting at the same space points 
for a smooth transfer of atoms. The adiabaticity in this 
case is limited by the smallest potential curvature, that 
is, the curvature of Vol (Eq. (3)) along the radial (y 
or z) direction, because the typical laser waist satisfies 
w > 4Aol- The separation between the ground state and 
the first excited state energies huj z gives a time scale 
2tt/uj z , where the angular frequency is defined as 

2 d 2 V L 



' ad 



dz 2 



z=0 



Wo 

„,,2 ' 



(15) 



Using the typical time scale tol = h/E r with the recoil 
energy E r = S 2 /c OL /2m of the optical lattice, the time 
scale is expressed as 



r (3) 
'ad 



V Vo 



10 



TOL 



which gives a time scale r 
typical parameter values (Appendix A. 2). 



(2) 
ad 



1-2 



(16) 

10 2 ~'Vs for 



5.3 Step 3 

Now, the atoms are trapped in the state-dependent 
optical lattice. According to our construction, the length 
between two nearest-neighbor apertures in the substrate 
is larger than a few times of the wavelength of laser 
beams. Then it is difficult to load two atoms in the near- 
est neighbor minima of the optical lattice, because the 
lattice constant is Aol/2. To collide the atoms, we must 
move the |0) (|1)) atom right (left) by n x A l/2, if the 
initial separation between the two atoms is 2n x Aol/2; 
n may be more than 5. 

The rough estimate of the execution time can be ob- 
tained by considering the adiabatic transport by one lat- 
tice constant Aol/2 along the x-axis in the usual optical 
lattice Eq. (3). The adiabaticity condition is then given 
by m(AoL/2To L ) 2 /2 ~ fiu x = hk oh ^2V /m. Thus, the 
adiabatic time scale for the motion with n lattice con- 
stant is estimated as 



t$ * n x Tol 



TOL, 



(17) 



where tol has been defined in Step 2. 

In our choice of the hyperfine states for |0) and |1), 
the state-dependent lattice potentials are given by Eq. 
(A-ll) in Appendix A. They are superpositions of two 
counterpropagating plane waves V± and there is an inter- 
ference between them. Then, the simple adiabatic argu- 
ment cannot apply to this situation. Actually, as shown 
in Appendix A, the fidelity is an oscillating function of 
the operation time for a given n; in other words, the 
adiabaticity time is not a monotonous function of n. We 
find through the numerical simulation that the execution 
time is roughly a few tens of the typical time scale tol ■ 

5.4 Step 4 

Let tpo( x ) be the ground state wave function of an 
atom in one of the optical lattice potential wells . If two 
atoms in the ground state are put in the same potential 
well, the interaction energy is given by 

2a s E r 



7rA 



OL 



ipQ(x)dx. 



(18) 



To acquire a phase tt necessary for the controlled- Z gate, 
for example, these atoms must be kept in the same po- 
tential well for £hoid = 7i7i/tAnt- 

For the s-wave scattering length a s — 5.19 nm between 
two 87 Rb atoms in \F = 1,itlf — 1) and \F — 2,rriF — 
l) 21 ) and typical parameters of the optical lattice (see 
Appendix A), we find Thoid of a few msec, which is the 
largest execution time among all the steps. This step does 
not involve any nonadiabatic process and we believe the 
operation is executed with a high fidelity by fine tuning 
the holding time £ n oid- 

5.5 Step 3', 2' and V 

It is clear that Step rC is an inverse process of Step 
n and time required for and the fidelity of the inverse 
process are the same as those for a forward process. The 
mathematical proof is described in Appendix B. 
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5. 6 Total execution time and fidelity 

We have estimated the time required for each step, un- 
der the assumption that each step gives a fidelity close to 
1. The overall execution time of a two-qubit gate is given 

by T ovoran ~ 2 (t^ } + t£> + r^) +r hold . The rough es- 
timate under the adiabaticity condition gives <~ 10 ms. 
Through the numerical simulation of the Schrodingcr 
equation with realistic parameters, as shown in Appendix 
A, we obtain 7.87 msec. 

The overall fidelity is estimated as follows. Each of 
Steps 1, 2, 3, 3', 2' and 1' involves two independent pro- 
cesses. For example, sending an atom from the vicinity of 
the substrate to the space point in Step 1 is realized with 
the fidelity 0.99 for each atom. Since this step involves 
two atoms, the fidelity associated with this step must be 
0.99 2 . Thus the overall fidelity is 0.99 12 - 0.886, where 
we assumed that the interaction time is tunable so that 
Step 4 gives a fidelity arbitrarily close to 1. The fidelity 
may be improved by spending more time for each step. 

Execution time may be shortened by applying stronger 
laser fields, because the separation between neighboring 
energy levels is enlarged by this, which makes fast adi- 
abatic operation possible. Besides, stronger laser fields 
suppress tunneling of an atom in the optical lattice to 
the adjacent potential minima. The ground state wave 
function is squeezed in a stronger lattice potential well, 
leading to a larger value of J ip 4 dx. Then, Step 4 is also 
executed in a shorter time. To make this statement more 
concrete, let us make a harmonic approximation to the 
potential well of the optical lattice. One can see that 
£hoid is reduced by a factor V 3 ^ 4 if Vq is made larger; 
for example, ihoid is multiplied by 0.178 if Vo becomes 
10 times larger. We expect that overall execution time 
can be reduced by one order of magnitude by employ- 
ing laser beams which are ten times stronger than the 
current value. 

In case the two atoms are in a general position, not 
necessarily along a primitive lattice vector of the optical 
lattice, we may introduce two orthogonal optical lattices 
as shown in Fig. 6 (a). Then the selected atomic states 
|0) of one atom and |1) of the other meet at the inter- 
section of the two optical lattices to acquire the extra 
dynamical phase. It should be noted that two-qubit gate 
operations may be applied simultaneously and indepen- 
dently on many pairs of qubits. 

6. Conclusion 

We proposed a scalable design of a neutral atom quan- 
tum computer with an on-demand interaction. A qubit 
is made of two hyperfine states of an atom. Associated 
with each atom, there is an optical fiber, through which a 
trap laser and a gate control laser are supplied. The lat- 
ter is used to implement one-qubit gates by two-photon 
Raman transitions. We have shown that a selective two- 
qubit gate operation is possible if atoms are trapped by 
an array of NFFD traps with a variable aperture size. 
This would allow for true two-qubit gates, rather than 
performing the same gate operation with many atoms in 
massive parallelism as in the experiments by Mandel et 
«Z. 9 ) We have obtained an uryner bound of the erate on- 



Fig. 6. (Color online) Two-qubit gate acts on atoms A and B 
in both cases. Two orthogonal optical lattices (light red contours) 
over two dimensional NFFD traps (dark black contours) . Atoms A 
and B are in arbitrary positions, not necessarily along one of the 
primitive vectors of the array. The screen is looked down from the 
above. 



eration time of 7.87 ms with the corresponding fidelity 
of 0.886. It may be possible to further reduce the ex- 
ecution time or increase the fidelity by taking different 
time-dependence of the control parameters. Reduction of 
the gate operation time is also possible by increasing the 
laser intensity. 

We believe that most of our proposal can be demon- 
strated within current technology. The severe technical 
challenge is to attach the ends of the optical fibers to 
the apertures of the substrate. We strongly believe that 
the rapid progress in nanotechnology makes it possible 
to solve this in the foreseeable future. Some other issues, 
such as trapping atoms near material surfaces, how to 
form an optical lattice in close proximity to a surface, and 
how to load individual traps with single atoms, should 
be addressed seriously toward realization of our proposal. 
We plan to study these issues in the future and would 
like to encourage experimentalists to design experiments 
based on our proposal. 
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Appendix A: Numerical Simulations of Two- 
Qubit Gate Operation 

The purpose of Appendix A is to estimate the execu- 
tion t.impt; nf ^t.prtt; 1 t.nrnncrb 3 in tho twrt-mihit. cffitp 
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operation introduced in Sec. 4.2 by solving the time- 
dependent Schrodinger equation 



ih— = -— v ip 

at 2m 



Vip 



(A-l) 



numerically. Using the experimentally relevant parame- 
ters, we calculate the optimized time to implement the 
gate operation keeping the fidelity 0.99 for each step. 
Here, the fidelity is defined as the overlap between the 
ground state wavefunction ip^ (x) at the final potential 
configuration in Steps 1, 2 and 3 and the solution of Eq. 
(A-l): 



dxrix,t)^{ X ). 



(A2) 



A.l Step 1 

A red-detuned laser beam passing through an aper- 
ture with a radius on the order of the laser wavelength 
forms an attractive potential U F {x,t) of Eq. (8) to an 
atom. The potential profile is a function of time if the 
aperture radius changes as a function of time. We solve 
Eq. (A-l) with V — U F (x,t) numerically to evaluate the 
time required to transfer an atom from the vicinity of an 
aperture to a space point, which will be a minimum of 
an optical lattice potential when it is turned on. 

We use the Di transition of 87 Rb atom to trap an atom 
and take the numerical values 

length Af 



energy 



time 



Tp = 

U 



795.118 nm 
k B x 20.1 yK 
h x 4.2 x 10 5 Hz 

0.38 //s, 



(A3) 



for the relevant parameters for definiteness, where the 
laser intensity is put to Iq = 10 5 W/cm 2 . 14 ) Here Af is 
the wavelength of the trap laser and Uq is the potential 
depth. The wavelength Af is fixed through the detuning 



ujl — ujq = 2nc 



1 

a7 



i 

An 



= -4.1 x 10 11 Hz, 



where Ao = 794.979 nm is the wavelength of the Di 
transition. 21 ' With these parameters, the lifetime of 
atoms due to the photon recoil heating is given as th ~ 
Uo/2Ep'j s ~ 1 sec with the recoil energy Ep = h 2 /2m\ F 
and the photon scattering rate 7 S = r e Sl 2 ff /4A 2 s , which 
is much longer than the total execution time. 

We rewrite the Schrodinger equation in a dimension- 
less form using the values given in (A-3), and intro- 
duce the cylindrical coordinates by taking advantage of 
the cylindrical symmetry around the z-axis. By writing 
<p = rip, the Schrodinger equation is put in the form 



d<P 



d 2 <p 

df 2 



ldcp 
f df 



+ U ¥ <P, (A-5) 



where dimensionless variables are denoted with a tilde 
and a = h 2 /2m\ 2 ? U ~ 2.2 x 10~ 4 . 

We take the initial aperture radius a; n i = 1.5 Af, cor- 
responding to z ini ~ 2.0, and the final aperture radius 
afi n = 2.8 Af, for which Zfi n ~ 8.0, in our computation 
Isee Fic\ 21. Let tW he the rhmensinnless time rermirerl 
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Fig. A-l. (Color online) Time evolution of the fidelity T for 
ft 1 ) = 1500 (a), 2000 (b), and 2500 (c). Since the fidelity increases 
monotonically with -fW, the execution time is evaluated at which 
the fidelity reaches 0.99 by interpolating the data, as shown in the 
inset. 



to change the aperture size from ai n i to ag n . The aperture 
size is changed as 

3(f) = a ini + (a fin - a ini ) sin 2 (^I)-*) (° < * < f 

(A-6) 

Here, we have vanishing initial and final velocities da/dt 
to avoid oscillatory behavior of the wave function during 
time evolution. 

We solved the Schrodinger equation (A-5) numerically 
and found = 2500, that is, tW = t^t f = 950 /its, 
is required to attain the fidelity 0.99 [see Fig. A-l]. This 
t^ 1 ) is consistent with the adiabaticity bound tAj in (14). 

A.2 Step 2 

We switch from the NFFD potential U F (x, t (1) ) to the 
optical lattice potential Uol(x) with the time scale 
as 



(A-4) V F (x,t) = cos 2 (^727*) u f(x,t {1) ) 



Voh(x,t) = sin 2 



2r( 2 ) 



(A-7) 



t) Uol(x) 



to this end. Here U F (x, r' 1 - 1 ) is the NFFD trap potential 
at t = and 

Uol(x) = -V cos 2 (2k OL x)e- 2( -y 2 +^- z ^ 2 ^ w2 (A-8) 

is the optical lattice potential, where z m is the z- 
coordinate of the minimum of the NFFD potential at 

It turned out to be convenient to employ different scal- 
ing for physical parameters here from those in Step 1. We 
use scales 



length A l 
energy E T 

time 



h 2 h 2 
2m 

h 



Tol = e: 



785 nm 
k B x 0.176 
hx3.7x 10 3 Hz 

43 /is 

(A-9) 
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before. 

Now the Schrodinger equation (AT) is put in a dimen- 
sionless form as 



.dip 



1 



dt 



(27T) 



; W 2 ^ + [V i J (x,t) + V F (x,t)}i> (A-10) 



For definiteness, we take V~o/E r — 40, Uo/E r = 114 and 
w/Xol = 4. We have taken Vo to be of the same order as 
that in Mandel et al. 9 ^ It should be noted that there is no 
rotational symmetry in the above Schrodinger equation 
and it must be solved by using the Cartesian coordinate 
system. 

We have solved the Schrodinger equation (A-10) nu- 
merically and found that the time scale = 560 /is is 
required to attain the fidelity 0.99 with the similar anal- 
ysis done in the Sec. A.l. Thus the above value of is 
consistent with this estimation. 

A. 3 Step 3 

The optical transitions shows that the effective poten- 
tials acting on |0) = \F = l,mp = 1) and |1) = \F = 
2, uiF = 1) are 



V w {x) = - A V + (x) + - 4 V-(x), 



(ATI) 



respectively. Here, we use the following decompositions 



|0) = 



3 1 
2' 2 



1 1 

2' 2 



V3 
2 



3 3 
2' 2 



1 

2'' 



V3 


3 1\ 


1 1\ 1 


3 3\ 


1 1 


2 


2' 2/ 


2' 2/ + 2 


2' 2/ 


2' 2 



(AT2) 

where 1 3/2 , 1/2) 1 1/2, 1/2) in the right hand side denotes 
a vector with nuclear spin (/ = 3/2, I z — 1/2) and elec- 
tron spin (S = 1/2, S z = 1/2), for example. If the po- 
larization angle 9 is increased, V\ ){x) (V|i)(x)) moves 
right (left). The choice of these hyperfine states gives an 
exactly counter-moving potential Eq. (ATI) for |0) and 
|1), while the states used by Mandel et al. 9 ^ yield the 
asymmetric potential given in Eq. (6). Hence we restrict 
ourselves to the motion of |1) state below. 

We take the number of the lattice constant over which 
an atom is transported to be n. In actual experiments, 
the angle may be changed by using an electro-optic 
modulator (EOM). The angle 9 is then limited within 
the range < 9 < it. This means that we have to reset 
9 to zero as soon as it reaches to 9 = n, if we need to 
change 6 beyond it (Fig. A-2). The trick is that even 
when 9 jumps by tt, cos# does not change at all. This 
reset must be done instantaneously so that the trapped 
atoms does not move in the optical lattice during the 
resetting time. By repeated use of this reset n — 1 times, 
it is possible to change 9 effectively from 9 — to 9 = nn. 



We have solved the Schrodinger equation 



chp_ 
dt 



2m 



vV + V|i)(t)v 



(A- 13) 




Fig. A-2. (a) Polarization must be rotated by 3tt to move an 
atom by three wavelengths of the optical lattice, (b) Actual change 
of 8 when EOM is in use. Even though 9 has discontinuity at t = Ti 
and T2, cos 6 is continuous in time. The effect of the angle 6 shown 
in (b) on an atom is the same as that of 8 in (a). 




Fig. A-3. (Color online) The fidelity as a function of t' 3 ) for 
n = 6. 



numerically. The angle 9 is changed as 
9(t) = nTrsin 2 ( 



(A-14) 



,2r( 3 ) 

in case an atom is transferred in the optical lattice by n 
site. We employ the same scaling as introduced in Step 
2, with which the Schordinger equation is put into a di- 
mensionlcss form as 



.dip 



1 




dt (2tt) 2 I dp 2 pdp ' p 2 ' d~z 2 J +V I 1 ^- 

(A- 15) 

Here, we use the cylindrical coordinates by taking ad- 
vantage of the cylindrical symmetry of the optical lattice 
around the z-axis. 

The Schrodinger equation is solved numerically and we 
find that the operation time , which gives the fidelity 
0.99, is 0.577 ms and 1.28 ms for two choices n = 3 and 
n = 6, respectively. Note that t^ 3 > is not a monotonous 
function of n. This is due to the fact that the lattice po- 
tentials V\q\ and V\i\ are superpositions of two counter- 



propagating plane waves V± and there is an interference 
between them. In fact, the fidelity is an oscillating func- 
tion of t*- 3 ' for a given n; an example is shown in Fig. 
A-3. These complicated properties make the execution 
time t^ 3 ) several times longer than the time scale 
obtained from adiabaticity requirement in Sec. 5.3. 
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Appendix B: Equivalence of JFo_>.t and Tt^o 

Let IV'o(O)) and \tpo(T)} be the eigenstatcs of the 
Hamiltonian at t = and t = T, respectively. Under 
the potential change V(t), and hence the Hamiltonian 
change H(t) 7 the state develops as 

|V(*)> = C/(t)|Vo(0)>, U(t)=Te-*ti H ^ dt . (Bl) 

The fidelity of this process is defined as 

.F„-t = (Vo(T)|V(T)) - (Vo(T)|C/(T)|Vo(0)). (B-2) 

For the reverse process H(T) — > H(0), we start with 
|^o(T")). The time development operator is = 
?7(t) _1 . Then the fidelity is 

Ft^o = (M0)\U(T)\MT)) 

= (Mt)\u(t)\mo))* 

= -^oUt = J^o^t, (B-3) 
where the reality of T is used. 
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